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Structured jets in TeV BL Lac objects and radiogalaxies. 

Implications on the observed properties 
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Abstract. TeV BL Lacertae objects require extreme relativistic bulk motion in the gamma-ray emission region, but at the 
VLBI scale their radio knots hardly move. The same sources show evidence, in radio, of a structure made of a fast spine plus a 
slow layer. We propose that this structure exists even on the spatial scale of regions responsible for the gamma-ray emission. 
One component sees the (beamed) radiation produced by the other, and this enhances the inverse Compton emission of both 
components. In addition, this allows the magnetic field to be nearly in equipartition with the emitting particles. The inverse 
Compton emission of the spine is anisotropic in its frame, possibly producing the deceleration of the spine by the Compton 
rocket effect. In this scenario, also the slow layer is a relatively strong high energy emitter, and thus radiogalaxies become 
potentially detectable by GLAST. 

Key words. Galaxies: jets — BL Lacertae objects: general — Radio continuum: galaxies — Radiation mechanisms: non- 
thermal — Gamma-rays: theory 



1. Introduction 

There is growing evidence from VLBI studies that pc scale 
jets in strong TeV BL Lacs move slowly (Edwards & Finer 
2002; Finer & Edwards 2004; Giroletti et al., 2004). On the 
contrary, the bright and rapidly variable TeV emission implies 
that at the jet scales where this emission originates, the jet 
should be highly relativistic. This is necessary in order to avoid 
the absorption of TeV photons by the IR radiation produced 
cospatially to the TeV emission (see e.g. Dondi & Ghisellini 
1995). Furthermore, fitting the SED of TeV sources with a sim- 
ple, one-zone homogeneous synchrotron self-Compton (SSC) 
model allows to determine a unique set of physical parame- 
ters of the emitting region (Tavecchio, Maraschi & Ghisellini 
1998), and in fact all authors applying SSC models derive sim- 
ilar values of the Doppler factor, in the range 10-20, when the 
TeV spectrum is not de-reddened by the absorption of TeV 
photons by the IR background (Tavecchio et al. 2001, Kino, 
Takahara & Kusunose 2002; Ghisellini, Celotti & Costamante 
2002; Katarzynski, Sol & Kus 2003), and larger (up to 50) 
when the TeV spectrum is dereddened (Krawczynski, Coppi 
& Aharonian 2002; Konopelko et al. 2003). It is therefore clear 
that the jet must suffer a severe deceleration from the 7-ray 
emitting zone ('^O.l pc from the jet apex) to the VLBI (~1 pc) 
scale. 
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Prompted by these observational evidences, Georgano- 
poulos & Kazanas (2003, hereafter GK03) have proposed that 
if the entire jet is rapidly decelerating in the 7-ray zone, then 
the base of the jet, still moving fast, will see the radiation 
produced at the end of the deceleration zone relativistically 
boosted. This "extra" radiation will favor the inverse Compton 
emission, allowing to derive less extreme values of the phys- 
ical parameters with respect to a pure one-zone SSC model. 
Our paper is germane to the one of GK03, but we study the 
alternative hypothesis that the jet is structured not in the radial 
direction, but in the transverse one, being composed by a slow 
layer and a fast spine. 

We are motivated by the recent observational evidence 
coming from detailed VLBI (including space VSOF observa- 
tions) radio maps, showing, in Mkn 501, a limb brightening 
morphology, interpreted as evidence of a slower external flow 
surrounding a faster spine (Giroletti et al. 2004). Similar results 
have been obtained for a few radiogalaxies (Swain et al. 1998; 
Owen et al. 1989; Giovannini et al. 1999). Apart from observa- 
tional evidence, a spine-layer configuration for the jet has been 
proposed in the past on the basis of theoretical arguments (e.g., 
Henri & Felletier 1991). In addition, the existence of a veloc- 
ity structure has also been suggested to explain some observed 
properties of radiogalaxies, such as their magnetic field config- 
uration (Komissarov 1990; Laing 1993), and to overcome prob- 
lems unifying radiogalaxies with BL Lac objects (Chiaberge et 
al. 2000). 
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We also hope that this assumption helps to find a possible 
reason for the deceleration of (at least a part of) the jet, which 
was postulated by GK03. A slow layer could in fact be the re- 
sult of the interaction of the "walls" of the jet with the ambient 
medium, or simply be the result of a jet acceleration which is 
a function of the angular distance from the jet axis, ' produc- 
ing a velocity structure. Kelvin-Helmoltz instabilities (for re- 
view see Ferrari 1998), while important for the formation of the 
layer, may not decelerate (and destroy) the entire jet (e.g. Bodo 
et al. 2003), especially if the jet itself is not continuous, but "in- 
termittent", as in the internal shock scenario (Ghisellini 1999, 
Spada et al. 2001, Guetta et al. 2004). Note in fact that the dy- 
namical timescale involved in internal shocks is approximately 
the light crossing time across the source, while instabilities can 
grow with the sound speed. 

As in the decelerating scenario proposed by GK03, there 
will be a strong radiative interplay and feedback between the 
layer and the spine: both parts see extra seed photons coming 
from the other part, and this will enhance the inverse Compton 
emission of both components. This may help explaining why 
also radiogalaxies can be relatively strong 7-ray emitters, as 
suggested by the Comptel and EGRET (onboard the Compton 
Gamma Ray Observatory) detection of Centaurus A (Steinle et 
al. 1998 and references therein), the recent identifications of 
NGC 6251 with an EGRET source (Mukherjee et al. 2002), 
and the possible detection of M87 at TeV energies (Aharonian 
et al. 2003). This emission, coming from the inner part of the 
relativistic jet of radiogalaxies, should be characterized by a 
pronounced variability: this could help to distinguish it from 
the high energy radiation coming from the more extended (kpc) 
parts of the jet, as suggested by Stawarz et al. (2003). 

In Section 2 we present the basic assumptions of the model 
and in Section 3 we discuss what we think is a major out- 
come of the scenario we are proposing, namely the possibil- 
ity that the spine can recoil under the effect of its own in- 
verse Compton emission (Compton rocket effect). We apply 
the model in Section 4 to Mkn 421, Mkn 501, Cen A and also 
to NGC 6251, assumed to be "paired" with the classic BL Lac 
PKS 0735H-178. Even if the presented "fits" ^ are not unique 
solutions (given the large number of free parameters), they il- 
lustrate the radiative feedback between the spine and the layer 
Since it is natural, in our model, that also the inverse Compton 
emission of the layer is enhanced, we stress in Section 5 that ra- 
diogalaxies, not only blazars, can be bright 7-ray emitters, and 
we present a preliminary list of possible candidates for detec- 
tion by GLAST In Section 6 we discuss our results and draw 
our conclusions. 



2. The model 

We assume that the layer can be approximated as a hollow 
cylinder, with external radius R2, internal radius R and width 



' Note that the possibility of structured jets has been proposed and 
somewhat explored also for Gamma-Ray Bursts, see Rossi, Lazzati & 
Rees 2001; Zhang & Meszaros 2002. 

^ We call "fit" what really is a comparison of the model against the 
observational data. 




Fig. 1. Cartoon illustrating the layer-nspine system. 



Ai?^ , as measured in the comoving frame of the layer ^. The 
comoving volume of the layer is then VJ = tt{R2 — i?^)Ai?;'. 
Also for the spine we assume a cylindrical geometry, with 
the same R and a width AR'^, as measured in the comoving 
frame of the spine. The active volume of the spine is then 
Vg — ttR^AR'^. Fig.^illustrates the assumed geometry. The 
Lorentz factors of the spine and of the layer are Tg and Ti, re- 
spectively, with c(3s and c/3; the corresponding velocities. Since 
the spine and the layer move with different Lorentz factors, the 
radiation emitted by the spine (layer) is seen boosted by the 
layer (spine). With respect to a comoving observer at the same 
distance from the spine (layer), the radiation energy density is 
enhanced by a factor ~ (T')^, with F' given by 



F' = TsTiil-PsPi) 



(1) 



Both structures emit by the synchrotron and the inverse 
Compton processes. The energy distribution of the emitting 
electrons, N{'-f), is assumed to extend down to 7niin, and is 
assumed to have the shape: 



iV(7) ^ K-f-"' 
Nil) - 0, 



lb 



e ^cut ; 7 > 7„i„ 

7 < 7min (2) 



The normalization (i.e. K) of this distribution is found by im- 
posing that A^(7) produces a given intrinsic synchrotron lumi- 
nosity, which is an input parameter of the model. 

The seed photons relevant for the scattering process are 
produced not only by the spine (layer) electrons, but also by the 
layer (spine) ones. There is a strong feedback between the two 
components, which determines the amount of inverse Compton 
radiation emitted by both structures (if this process is important 
with respect to the synchrotron one). As a general rule, this 
feedback increases the inverse Compton flux, since both the 
spine and the layer see an enhanced radiation energy density. 
Since the ratio between the radiation and the magnetic energy 
densities [/e/C^rad measures (even if it is not equal, for scat- 
terings in the Klein Nishina regime) the ratio between the in- 
verse Compton and the synchrotron luminosities, an enhanced 
f/iad in turn implies a larger magnetic field. Therefore both 
the synchrotron and the inverse Compton luminosities can be 
produced by a reduced number of relativistic electrons (radia- 
tively cooling in a shorter time with respect to the case of no 



Primed quantities are in the rest frame of the spine, double primed 
quantities are in the frame of the layer. 
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feedback). This fact bears important consequences for the en- 
ergetics and the dynamics of the jet, and it will be discussed in 
Section 3, while in Section 4 we give an illustrative example 
for the case of Mkn 421 . 

The length of the layer, as observed in a frame comoving 
with the spine, is R[ — li^ /V . Analogously, the length of 
the spine, as observed in the frame comoving with the layer, 
is i?g = R'g/T' . In the following, we will always assume that 
the layer is longer than the spine, even in the frame of the spine. 
To calculate the radiation energy density of one component as 
observed by the other, we will assume the following: 

- Let us call R = {R2 + i?)/2. In the comoving frame of the 
layer, the radiation energy density in the entire cylinder is 
assumed to be U'l = L\ /{ttR?c). In the frame of the spine, 
this radiation energy density is assumed to be boosted by a 
factor (r')^ i.e. U'l = {T'fU'i. 

- In the comoving frame of the spine, the radiation energy 
density within the spine is assumed to be U'^ = L'^ / {n R'^ c) . 
In the frame of the layer, this radiation energy density is 
observed to be boosted by (F')^, but also diluted (since the 
layer is larger than the spine) by the factor Ai?g/Ai?; — 
{AR'JT')/Ali;. 

While the synchrotron and the SSC emission in the comov- 
ing frame is assumed to be isotropic, the inverse Compton pro- 
cess between electrons of the layer (spine) and seed photons 
produced by the spine (layer) is highly anisotropic. Dermer 
( 1 995) found the pattern of the emitted radiation from a moving 
blob immersed in a bath of seed photons (e.g. corresponding to 
the radiation produced by the broad line region of a powerful 
blazar), and pointed out the fact that in this case the external 
Compton radiation is more beamed than the synchrotron and 
SSC emission. 

In our case the component contributing to the "external" 
radiation is not at rest with the distant observer, but moves. To 
find out the pattern of the emitted radiation it is convenient to 
move to the comoving frame of the emitter of the seed pho- 
tons. Consider then the seed photons produced by the layer, 
and an observer comoving with the layer. In this frame, the 
spine is moving with F', and the photon frequencies produced 
by the spine are blueshifted by the Doppler factor * 6s.i- Going 
to the frame of the distant observer, these photons are further 
blushifted by (5; . But the distant observer will see the same pho- 
tons blushifted by Ss- This implies 



Ss,i Si 



(3) 



We can repeat the same argument for photons produced by the 
layer. We then have 



Si 



1 



(4) 



This nice argument is due to GK03. 

In the frame of the layer the external Compton radiation 



produced by the spine follows a pattern oc (5j^|^" (Dermer 

"* (5 = r^^(l — f3 cos 9)^^, where 9 is the viewing angle. Note that 
r and 9 are different in different frames. S^^i is defined as the beaming 
factor of the radiation produced in the spine as observed in the layer. 



1995), where a is the spectral index of the emission [F(i/) oc 
i^^"], while the synchrotron and SSC emission follows the 
usual pattern oc (5^^". If I'{t^') is the monochromatic intrin- 
sic intensity produced by the spine, we have 



/(^) = /'(^')<r"'5f+" = l'ii^')S'+- (^1) "(EC) 
I{i^) = /'(i^')'5f+" (S, SSC) (5) 

Here S stands for synchrotron, EC for "external Compton" 
(scattering of the seed photons coming from the layer). In the 
case of radiation produced by the layer, the transformation is 
the same, with 5s and 5i interchanged. 

Free parameters — For each component the input parameters 

are: AR', B, L|„j, F, 7min, 7b, 7cut, ^1,712. Parameters equal 
for both components are R and the viewing angle 6. The outer 
radius of the layer i?2 must also be specified, but we will always 
use i?2 ~ 1.2R. For 71,2 > 3, 7cut becomes unimportant and 
we have a total of 18 parameters. While they appear (and are) 
many, we stress that the aim of this paper is to discuss the main 
effects of having the layer-nspine structure, and not (yet) the ex- 
act determination of the physical quantities inside the source. 
In other words, we have to model two structures by observing 
the radiation coming from only one of them, and this will leave 
some ambiguity, unless we find information about the layer 
(spine) even if we are observing the spine (layer) emission ^. 

3. Energetics and dynamics 

As we will see in the next Section, the power spent by the jet to 
emit the inverse Compton radiation is significant, if compared 
with the total bulk kinetic power contained in the relativistic 
electrons and protons. As described above, this luminosity is 
emitted anisotropically in the comoving frame of the spine, 
hence the spine must recoil. Since this Compton rocket effect 
is a potentially important deceleration mechanism, we discuss 
it here in more detail. 

The most convenient frame where to study the dynamics of 
the spine is the frame comoving with the layer. In this frame the 
synchrotron radiation produced by the layer is isotropic, with 
energy density t/" „. To avoid an excess of symbolism, from 
now on (unless otherwise noted) the unprimed quantities are 
measured in the frame of the layer, and the primed quantities 
are evaluated in the frame of the spine. For the same reason we 
let the random Lorentz factor unprimed, with the notion that 
these Lorentz factors are measured in the frame of the spine, 
where they are assumed to be isotropically distributed. 

The total Lorentz factor (7) of the electrons is the superpo- 
sition of the relative bulk (F) and random (7) Lorentz factors. 
Assume that the random velocity forms an angle 9' with respect 
to the jet axis, in the spine frame. From Rybicki & Lightman 



One possibiUty could be the observation, in radiogalaxies, of 
emission lines resulting from photoionization due to the spine radia- 
tion (Morganti et al. 1992), or the knowledge of the total power carried 
by the entire jet through estimates of the total energy and age of the 
extended radio structure (as done by Rawlings & Saunders 1991). 
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(1979) we have: 

^ P' cos e' + (3buik 
~ 1 + Pbuikl3'cos0' 
13' sin 0' 
"" " rbuik{l + (3buik(3'cos0') 

The total 7-^ is 

= (l + /36„ifc/3'cos0')Vr' 

If the particle distribution is isotropic in the 
average over angles gives 

^ J2nsm0'j^{0')d0' 

^^ ' 47r 

= ^ J {l+PbuikP'l^'fdii' 

, iPbulkl^') 

3 



(6) 



(7) 
spine frame, the 



^2r2 



(8) 



the factor in the square parenthesis becomes (4/3) for ultrarel- 
ativistic speeds. 

The loss of energy of the jet is proportional to (7^), and the 
loss of momentum is the component along the jet axis of the 
loss of energy. In other words, we have to calculate the quantity 



^2^^,'^\0')d^^' 2 

(7 )z = -^ = -^PbulkPl 1 



(9) 



Now assume that the spine carries Np protons (total num- 
ber) and N(. leptons. The loss of momentum of the jet is de- 
scribed by: 



dr_ ^ 4 aTcNeUsyn{r)z 

dt 3 NpnipC^ + Ni,{'^)meC^ 

8 aTcA^et/syn(7')r' 

9 NpfUpC^ + Ne{j)mec'^ 



We remind that the bulk Lorentz in Eq. dlOt is the relative 
Lorentz factor (which is the Lorentz factor of the spine as mea- 
sured in the layer). A detailed analysis of this Compton rocket 
effect in AGN jets has been presented by Sikora et al. (1996). 
They derive the equations regulating the drag effect and ap- 
ply the results to the case of jets in radio-loud quasars with 
strong emission lines. Eq. (llOt . although obtained in a rather 
simplified way, coincides with the results of Sikora et al. (for 
the specific case of isotropic radiation field). 

To find the value of Tg as measured by an observer at earth, 
one must invert Eq. ^ It must be noted that, in general, (7} 
and (7^) are not constant, but can change due to e.g. radiative 
cooling. We have integrated Eg. 1101 assuming: 

- The particle distribution remains unaltered for a time equal 
to the light crossing time of the spine (i.e. ARs/c) as mea- 
sured in the frame of the layer. This corresponds to assume 
that the time of the injection of particles throughout the 
spine lasts for a similar time. 
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Fig. 2. Example of the SED produced by the spine-layer sys- 
tem, using the parameters listed in Tab. 1 . Dashed lines corre- 
spond to the emission of the spine (layer) without taking into 
account the seed photons coming from the layer (spine). Data 
from Plan et al. (1998) and Djannati-Atai et al. (1999). 



- After this time we calculate the actual 7cooi at the given 
time, and assume that the particle distribution vanishes for 
7 > 7cooi- We neglect adiabatic losses, for simplicity. 

- We then calculate the new values of (7) and (7^) and the 
new value of F. 

- We invert Eq[l]and finally find Vs. 



4. Results 

^10) 4. r. Fimng Mkn 501, Mkn 421 and Cen A 



We here apply the model to the two best studied TeV BL Lacs: 
Mkn 501 and Mkn 421, and to the radiogalaxy Cen A. 

We stress that our aim is not to precisely fit the SED of 
these sources, since our results are not unique. For instance, 
the parameters for the layer of Mkn 421 and Mkn 501 have 
been chosen with the aim to increase the allowed value of the 
magnetic field of the coiTesponding spines. In other words, we 
applied a theoretical prejudice, aiming to bring these sources 
closer to equipartition than a simple one zone SSC model al- 
lows to. Other poorly constrained parameters are the width of 
the layer and its bulk Lorentz factor. For the former we simply 
require AR" > R, for the latter we have the requirement that 
the counterjet is, in these sources, invisible up to the few tens 
of m.a.s. scale, leading to a limit of F; > 2-3. 

In conclusion, our aim is to demonstrate that our model can 
consistently work both for radiogalaxies and BL Lacs, and this 
allows to draw some general conclusions about their jets and 
the importance of the inverse Compton emission. The input pa- 
rameters used for the models shown in Figs. 2-5 are all listed in 
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Fig. 3. The SED of Cen A is modelled by the spine-layer sys- 
tem, using the parameters listed in Tab. 1 . Dashed lines corre- 
spond to the emission of the spine (layer) without taking into 
account the seed photons coming from the layer (spine). For 
the spine the dashed and continuous lines overlap. Data from 
Chiaberge, Capetti & Celotti (2001) and references therein. 



Fig. 4. Example of the SED produced by the spine-layer sys- 
tem, using the parameters listed in Tab. 1 . Dashed lines corre- 
spond to the emission of the spine (layer) without taking into 
account the seed photons coming from the layer (spine). For the 
data points see Costamante & Ghisellini (2002) and references 
therein. 



Table 1 , while Table 2 reports the values of the kinetic powers 
carried by the spine and the layer, the average (and the aver- 
age of the square) of the random Lorentz factor of the emitting 
electrons, and the relative bulk Lorentz factor F'. The kinetic 
powers are defined as 



U 



ttR^T^PcU^ 



(11) 



where the energy density Ui refers to: electrons (C/g = 
J N{'y)'ymeC^djy, magnetic field (Ub — B^/Stt); protons 
(Up), assuming one cold proton per electron; total radiation 
(C/rad) and synchrotron radiation ([/gyn)- 

In all cases shown in Figs. 2-5, continuous lines show the 
results of the models including the spine-layer radiative feed- 
back, while the dashed lines are the results of the same models 
neglecting the feedback. 

- Mkn 501: Fig. 2 shows the "fits" to the SED of Mkn 501 
in its flaring state of April 1997 (Plan et al. 1998). Note 
the rather moderate bulk Lorentz factor used for the spine 
(Ts = 15), which has the same value of the one used by 
GK03. The rather large value of the magnetic field ensures 
nearly equipartition between magnetic and particle energy. 

- Mkn 421: Fig. 3 shows our model compared with different 
simultaneous SED of this source, and aiming to "fit" the 
highest state. The bulk Lorentz factor of the spine is F^ — 
17, which at a viewing angle 6 — 2.5° leads to a beaming 
factor 6 — 22. The remarkable difference with respect to 
a simple one zone SSC models is the much larger value of 
the magnetic field: 1 Gauss compared with B ^ 0.02-0.1 
Gauss (see the references quoted in the Introduction). This 



brings the magnetic energy in equipartition with the par- 
ticle energy, dominated in this case by the energetic elec- 
trons, which have a very large average random Lorentz fac- 
tor ((7) = 1900, and therefore {j)me ~ rup). This is due 
both to the relatively large value of 7niin and the flat elec- 
tron slope at low energies, in turn required in order not to 
overproduce the IR emission. 
- Cen A: As can be seen in Fig.|3] the synchrotron spectral 
component in this object is particularly naiTow, with a peak 
in the far IR. This indicates a coiTesponding narrow energy 
distribution of the emitting electrons. We have assumed 
a viewing angle of 6' = 40°, corresponding to a modest 
beaming of the layer emission (S = 1.2), and a severe de- 
beaming of the spine radiation (S = 0.28). Comparing the 
jet powers of this source with Mkn 421 or Mkn 501, we 
find that the jet of Cen A is more powerful, with most of 
the power carried by the spine. 

4.2. "Pairing" BL Lacs witti radiogalaxies 

We here try to see if our spine-nlayer structure can explain at 
the same time the emission from a "classical" BL Lac object, 
PKS 0735H-178, and the radiation observed for the radiogalaxy 
NGC 625 1 , which has been recently associated with an EGRET 
source. In other words, we check if these two apparently very 
different objects, at very different distances (but both emitting 
high energy 7-rays), can be considered as "paired". Note that 
the radio luminosity Li 4 at 1.4 GHz of the two objects, indica- 
tive of the extended (and unbeamed) radio power, is almost 
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Fig. 5. Example of the SED produced by the spine-layer sys- 
tem, using the parameters listed in Tab. 1 . The viewing angle is 
6 = 2° for the example shown in the upper panel, and 9 = 22° 
for the bottom panel. The model is compared with the data 
of the BL Lac PKS 0735+178 (top) and the FR 1 radiogalaxy 
NGC 6251 (bottom). Dashed lines correspond to the emission 
of the spine (layer) without taking into account the seed pho- 
tons coming from the layer (spine). Data from Ghisellini et al. 
(1998) and Chiaberge et al. (2003) and references therein. 

equal (for PKS 0735+178: L1.4 = 9 x lO^i erg s'^ Hz"!, 
Cassaro et al. 1999; for NGC 6251: L1.4 = 3 x lO^i erg s"! 
Hz~^, Laing, Riley & Longair, 1983, extrapolating from 178 
MHz with a radio spectral index a = 0.7). In this case the dif- 
ferences in the observed nuclear SED are due entirely to a dif- 
ferent viewing angle 9, enhancing the spine or the layer emis- 
sion for small and large 9, respectively. As a consequence, we 
have now much less freedom than before for the fitting, but, of 
course, now it is the choice of these two sources which is some- 
what arbitrary, even if both are 7-ray emitters and presumably 
both belong to the FR I class of radiogalaxies. 

- PKS 0735+178: This "classic" BL Lac has been detected 
by EGRET (Hartman et al. 1999), and shows a relatively 
flat X-ray spectrum, signature of the fact that at these en- 
ergies the inverse Compton emission is already dominat- 
ing. The synchrotron peak frequency lies in the IR-optical 
and possibly "moves" within these bands. The data shown 
in Fig. 4 are not simultaneous, and therefore give a rough 
idea of the entire SED. Nevertheless it is clear that this 
source belongs to the LBL (Low energy peak BL Lac, see 
Giommi & Padovani 1995) category. There is only a lower 
limit for its redshift, (z > 0.424); for the modeling we as- 
sume z = 0.424. The implied energetics are two orders of 
magnitude larger than for the previous sources, if a sim- 
ilar value of the Doppler beaming factor is appropriate. 



Accordingly, we have assumed a larger intrinsic power, and 
a larger value of the magnetic field. This brings the mag- 
netic energy closer to equipartition (with respect to a pure 
SSC model) with the energy contained in the relativistic 
electrons. 

NGC 6251: The layer assumed to be in the jet of PKS 
0735+178, if observed at 20°, becomes the dominant con- 
tributor to the SED and give rise to the spectrum of the 
nearby radiogalaxy NGC 6251 (z = 0.0249), as illustrated 
in the bottom panel of Fig. 4. Note that also in this case 
the data points are not simultaneous, and can give only a 
rough idea of the overall (nuclear) SED of this object. The 
beamed non-thermal component, being variable with a rel- 
atively large amplitude, could dominate only occasionaly 
the total flux, letting in other cases the component produced 
by the accretion disk be dominant. This could also explain 
why, in some observations, a relatively strong fluorescent 
broad iron line (at 6.4 keV) is visible (Gliozzi et al. 2004). 
Bearing in mind these uncertainties in the relative strength 
of the beamed vs accretion components, the slight overpro- 
duction of X-rays predicted by the model is not (yet) neces- 
sarily a failure of the model. Note also that when the layer is 
not illuminated by the spine (and this can happen, since the 
spine can be constituted by discontinuous blobs), then the 
predicted emission coiTesponds to the dashed line, which is 
closer to the X-ray data, but underpredicts the 7-ray flux. 
This may originate two different states of the source, and 
EGRET could detect the object only when the layer was il- 
luminated by the spine. Note that this source showed low 
amplitude, relatively fast ('^ 10** s) variability in hard X- 
rays (above 0.8 keV), in agreement with the idea that the 
non-thermal beamed continuum rivals the contribution pro- 
duced by the accretion disk, with the latter better visible in 
the softer X-ray band (Gliozzi et al. 2004). 



4.3. Jet deceleration by the Compton rocket effect 



As anticipated above, from Table 2 we can infer that the power 
released as inverse Compton radiation by the spines of Mkn 
501 and Mkn 421 is comparable to the total kinetic power car- 
ried by the jet. In this circumstances the Compton rocket effect 
is important to determine the dynamics. For illustration, Fig.|6l 
shows the evolution of the bulk Lorentz factor in the case of 
Mkn 421, calculated along the lines discussed in Section 3, and 
assuming the same parameters used for the fitting described 
above and listed in Tab. 1. We show the bulk Lorentz factor 
as measured by the layer and also by the observer at earth, but 
both as a function of the distance as measured in the frame 
of the layer We also show the logarithm of 7cooi- As can be 
seen, the recoil is very significant, decelerating the spine from 
Ts = 17 to Ts = 7.5 wifliin flie layer. Note that L = 3 is 
the bulk Lorentz of the layer, which is therefore the minimum 
value which the spine could attain. 
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Table 1. Input parameters of the models for the layer and the spine shown in Figg. |2 
viewing angles for PKS 0735+178 and for NCG 6251 are 3° and 20° respectively. 
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Table 2. Derived jet powers and parameters of the models for the layer and the spine shown in Figg.|2]-|5] 
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Fig. 6. The bulk Lorentz factor as a function of distance as mea- 
sured in the frame of the layer and in the earth frame. Also 
shown, as labelled, is the logarithm of 7cooi- We have adopted 
the same parameters as the ones used for Mkn 42 1 and reported 
in Table 1 . 



4.4. General outcomes and implications 

From the exercise of "fitting" the SED of the previous few BL 
Lacs and radiogalaxies, we can draw the following general con- 
clusions; 



- The feedback between the layer and the spine enhances 
the inverse Compton emission in both components. In the 
shown examples, the enhancement is around one order of 
magnitude, or even greater. 

- The magnetic field is larger, in this spine/layer scenario, and 
is consistent with being in equipartition with the relativistic 
electrons. 

- The larger magnetic field implies that, to produce the same 
amount of synchrotron radiation, less electrons are needed. 
If there is one proton per electrons, we need also less pro- 
tons. This implies a decrease in the total jet kinetic power 
(with respect to a pure SSC model). 

- The Compton rocket effect can influence the dynamics of 
the spine, since the external Compton luminosity (being 
anisotropic in the comoving frame) becomes comparable 
with the total kinetic power carried by the spine itself. 

- The examples shown here assume that the spine is active 
only when inside the layer. But other possibilities exist: for 
instance, the layer could be particularly narrow (i.e. small 
ARi), and allow the spine to survive the passage through 
itself (i.e. the electrons may not cool completely). In such 
a case, the spine will continue to emit the same amount of 
synchrotron radiation (if the particle distribution does not 
change), while the amount of high energy emission will 
drastically decrease. This might explain the "orphan TeV 
flares" (i.e. flares in the TeV band not accompanied by si- 
multaneous flares in the X-ray band) observed in simulta- 
neous RXTE/TeV observations (Krawczynski et al. 2004). 

- To be effective, the layer must of course be located in the 
same region of the jet where the spine undergoes dissipa- 
tion and emits. In other words, it must be compact, even if 
its width can be larger (factor ~10) than the spine width 
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(i.e. AR'i ^ lOAi?^). This implies variability of the layer. 
Even if its particle distribution is steady, in fact, its high en- 
ergy emission will change as the spine is illuminating it or 
not (e.g. compare the continuous and dashed lines for the 
layer in Figg. 2-5). 



5. GeV and TeV radiogalaxies? 

Also the slow layer produces a lot of GeV radiation, which 
remains visible even at large viewing angles (consider that 
9 — l/Fiayer ~ 20°). Therefore also radiogalaxies should be 
high energy emitters. Indeed, the three radiogalaxies detected 
at high energies so far, shown in Fig. show the characteris- 
tic double bump SED typical of blazars, suggesting a similar 
origin of their emission. This is more clear in Fig. |8l which 
compares directly the averaged SED of blazars with the SED 
of Cen A, NGC 625 1 and M 87 . 

For Cen A and NGC 625 1 , which have been identified with 
EGRET surces, we measure a ratio between the EGRET (lyF^) 
flux at 100 MeV and the radio flux of the core at 5 GHz which is 
60 and 300, respectively. This large ratio is also consistent with 
the TeV detection of M87, if the peak of its 7-ray emission lies 
between the EGRET and the TeV band. 

If all FR I radiogalaxies have a similar 7-ray to radio flux 
ratio, we can identify the best candidates for detection in fu- 
ture AGILE and GLAST observations. To this aim, we have 
averaged the ratio of the 7-ray to radio flux for the three radio- 
galaxies detected in 7-rays which results in 



logi^^i^^ = (2±0.5)logi.RFR, 



(12) 



where i^R.coic is the radio flux of the core at 5 GHz (Giovannini 
et al. 1998; Chiaberge, Capetti & Celotti 1999 and references 
therein), and F-y is the 7-ray flux at 100 MeV. Using Eq.ll2las 
an empirical guide to predict the 7-ray fluxes of radiogalaxies, 
we list in Tab. 3 the FR I radiogalaxies belonging to the 3C 
sample sample, in order of decreasing nuclear radio power, and 
the coiTesponding predicted 100 MeV flux. Assuming a sensi- 
tivity limit of GLAST of ^^ 5 x 10~^^ erg cm^^ s^^ above 
100 MeV for one year of exposure time, i.e. a factor 20 better 
than EGRET, (see e.g.: http://glast.gsfc.nasa.govi, we have that 
more than a dozen FR I radiogalaxies can be detected. 



6. Discussion 

We have explored the radiative and dynamical consequences 
of a structured jet, in which a slow jet layer is cospatial to a 
fast jet spine. The motivations for this study are primarily ob- 
servational, since radio data have recently found slow or null 
proper motion for the parces scale radio knots in TeV emitting 
BL Lacs. In addition, detailed VLBI radio maps show hints of 
a limb brightening for the jet in Mkn 501 . 

The jet may be born as structured, or it may be born with 
equal layer/spine velocities, with the layer being decelerated in 
the first ~ 10^^ cm from the central black hole. As long as the 
layer is dissipative, our results are independent of the mecha- 
nism producing a spine-nlayer structure. However, we note that 



Name 


log FR„coro 


log u^Fj 




(5 GHz) 
erg cm~^ s~^ Hz^^ 


(100 MeV) 
erg cm^^ s^^ 


3C84 


-21.37 


- 9.67±0.5 


3C274 
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-10.70±0.5 
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Table 3. Radio core fluxes at 5 GHz for 3CR FR I and pre- 
dicted gamma-ray fluxes. See the text for the assumptions. 



the accretion disk in low power BL Lacs could be character- 
ized by a low-efficiency accretion mode, and one of the pro- 
posed accretion disk solution in this case is the ADIOS struc- 
ture (Blandford & Begelman 1999), predicting that a sizeable 
fraction of the accreting mass at large radii, instead to inward 
spiralling, leaves the disk in the form of an outflow. In this case, 
in the close vicinity of the black hole, the ambient medium 
should be relatively dense, favoring the interaction of the jet 
walls with the external medium. This may cause the formation 
of the slow layer, and at the same time this interaction may be 
the primary cause for the transformation of the kinetic energy 
of the layer into random energy and then radiation. 

A low radiative efficiency accretion disk could also explain 
why the broad line emission is small in BL Lac objects and in 
FR I radiogalaxies. Further evidence of a change in the accre- 
tion mode between FR I and FR II radiogalaxies is provided 
by the "dividing line" between these two classes of sources in 
the radio-luminosity-host galaxy optical luminosity plane, as 
discussed by GhiselHni & Celotti (2001). 

If indeed the jet forms (or is born with) a layer-nspine struc- 
ture, then there is a radiative interplay or feedback between 
the two parts: each component would see an enhanced radia- 
tion field coming from the other component. This would in- 
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Fig. 7. The SED of NGC 6251, Cen A and M87. The data of 
NGC 625 1 and M87, for clarity, have been vertically shifted by 
the labelled amount. Note that the SED of these radiogalaxies 
show the same two-bump structure of blazars, as illustrated 
also in Fig.|8] Data for M87 are taken from Reimer, Protheroe 
& Donea (2004) and reference therein, except for the X-rays 
(Marshall et al. 2002), and the EGRET upper limit (Fichtel et 
al. 1994). 



evitably boost the inverse Compton radiation with respect to a 
completely homogeneous jet. 

It is also quite clear that, for an observed Compton to syn- 
chrotron power ratio, the fact that the radiation field is en- 
hanced also implies an increase of the magnetic field, with re- 
spect to an homogeneous source. This can solve an otherwise 
puzzling characteristic of HBL in general and TeV BL Lacs in 
particular: fitted with an homogeneous SSC model, they turn 
out to have very small magnetic fields, which are more under- 
equipartition (with the emitting particles) than in other blazars. 
Here we can fit the observed spectra equally well (not surpris- 
ingly, given that the free parameters are more than in the ho- 
mogeneous SSC model) with equipartition magnetic fields. 

Another important consequence of having greater magnetic 
fields is that a smaller number of electrons can produce the 
observed SED. This means that the global energetic demand of 
the jet is reduced with respect to a homogeneous SSC model. 

The fact that the dominant inverse Compton radiation is 
through scattering with "external" photons implies that the 
emission is highly anisotropic also in the comoving frame of 
the spine. To conserve momentum, the emitting spine must 
recoil and therefore decelerate. This is a manifestation of the 
"Compton rocket" effect, studied in the early eighties (see 
e.g. O'Dell 1981) as the means to radiatively accelerate jets. 
Somewhat ironically, we have shown here that this process can 
be important for the opposite reason. There is a precise link 
between TeV emission and deceleration, due to two reasons: 



Fig. 8. The SEDs of NGC 6251, Cen A and M87 are compared 
with the blazar sequence, as proposed by Fossati et al. (1998). 
The hard X-ray [2-10 keV] spectra of blazars come from the 
work of Donato et al. (2001). 



i) TeV emitting BL Lacs have the least powerful jets, and yet 
they move with bulk Lorentz factor equal or greater than the 
ones of other blazars; ii) to produce a significant TeV radia- 
tion, the mean energy of the emitting electrons must be large. 
Then in these sources we have (7)me ~ nip: the power caiTied 
by the jet in the form of protons and electrons is similar. We 
then propose that the jet deceleration at small (sub-pc) scales 
is more efficient in low power jets emitting high energy radia- 
tion. Consider also that the initially fastest sources are the ones 
suffering the most severe Compton rocket effect. 

Also the inverse Compton emission from the layer is en- 
hanced by the extra seed photons coming from the spine. This 
could be the reason why also radiogalaxies are relatively strong 
7-ray emitters. If this is the reason, then the layer and the spine 
must be cospatial, and therefore the 7-ray emitting layer must 
be located at ^100 Schwarzchild radii, as the spine. The 7- 
ray flux observed in radiogalaxies should then be variable, with 
timesc ales of the order of ivar ~ {R/c)/Si ~ a day or less, very 
similar (albeit somewhat longer, due to the smaller Doppler 
factor) to the typical variability timescale in blazars. We pro- 
pose this as a crucial test for our scenario. Note that Cen A is 
already known to vary with short timescales (0.5^ days) in the 
7-ray band (Kinzer et al. 1995; Steinle et al. 1998). 

A few radiogalaxies have been akeady detected at high en- 
ergies. The SED of their nuclear emission shows the charac- 
teristic double peak characteristic of blazars, and like these 
sources they can be "Compton dominated" (namely, the high 
energy component is more luminous than the synchrotron com- 
ponent). These general features are quite easily explained in the 
layer-Hspine scenario. We have then tried to predict which are 
the best candidate radiogalaxies to be detected by the future 
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GLAST mission, assuming a twenty-fold increase in sensitiv- 
ity with respect to EGRET. As a result, more than a dozen PR 1 
radiogalaxies should be detectable by GLAST, if their radio to 
100 MeV flux ratio is similar to the one of the three objects al- 
ready identified by EGRET and by TeV Cherenkov telescopes. 
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